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Abstract 

We investigate the spin accumulation effect in eccentric semiconductor multichannel rings with Rashba 
spin-orbit interaction and threaded by a magnetic flux. Due to the finite eccentricity, the spin polarization 
induced at the borders of the sample is anisotropic and exhibits different patterns and intensities at specific 
angular directions. This effect, reminiscent of the spin polarization drift induced by the application of an 
in plane electric field, could be used to manipulate and functionalize the spin polarization in electronic 
nanorings. 
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I. INTRODUCTION 



The spin manipulation and detection in all electrical devices have been the subjects of an im- 
pressive amount of research in the last years 1 . In particular semiconductors in which the spin orbit 
(SO) interaction plays a significant role, are the main candidates to accomplish such challenges. 
Among them, ring shaped geometries are particularly attractive to analize different interference 
phenomena in the presence of SO effects. As an example, Aharonov Bohm (AB) rings with uni- 
form SO interaction have been proposed as spin interference devices^ to explore the non trivial 
spin dependent Aharonov-Casher (AC) phase^ 4 -^. Indeed, in recent years the AC effect has been 
measured in a series of transport experiments on AB conductance oscillations performed in semi- 
conductor rings, for different intensities of the SO interaction^^. 

So far most of the theoretical analysis of SO effects in quantum rings have been restricted to ID 
geometrie o 10 ' 11112 ' 13114 - 15 . However realistic rings have a multichannel nature and many interesting 
phenomena rely on this fact. As an example, in Ref [16| it was shown that in a AB multichannel 
ring with Rashba SO coupling— a spin accumulation effect develops at the borders of the sample. 
Eventhough for an even number of electrons, a finite spin polarization in the direction perpen- 
dicular to the plane of the ring is generated and can be controlled with the magnetic flux. This 
phenomena, although sharing some analogies with the intrinsic Spin Hall Effect (SHE) studied 
in bar or strip geometries designed on 2DE G 18 ' 19i2 ° , does not require neither external currents nor 
electric fields or voltage drops applied 2 ^ 2 .. 

Actually in real semiconductor quantum rings (QRs) imperfections of the structure often oc- 
curs. Micrograph views of GaAs QRs suggests that real rings may present some degree of 
eccentricity 2 ^. In addition, recent advances in oxidation lithography enables the fabrication of 
semiconductor QRs 24 in which non perfect symmetric structures can be constructed. 

The effect of a finite eccentricity has been considered in previous studies of energy spectrum 
and electric polarization in QRs 25 26 . However these works have not taken into account the SO 
interaction, which is particularly strong in most of the QRs heterostructures. 

The goal of the present work is to show that by a combined effect of the Rashba SO interaction 
and the finite eccentricity, the spin accumulation develops anisotropic patterns along the sample. 
This effect, that shall be described in detail later on, will lead to different intensities of the total 
spin polarization in the angular direction. The typical spatial separation of the anisotropic patterns 
may render possible to sense the effect employing usual magneto optical detection techniques 27 . 
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The article is organized as follows. Sec|II] introduces the model hamiltonian for the eccentric 
QR with Rashba SO interaction. After giving some details of the calculations, the energy spectrum 
and eigenf unctions are obtained for different magnetic fluxes. These results are analyzed and 
compared to the ones of the concentric QR. In SecQlI]the spin polarization is computed and the 
different anisotropic patterns in the spin accumulation effect are analyzed in detailed. The last 
section |IV] is devoted to the conclusions and to elaborate on the experimental feasibility to detect 
the anisotropic profiles in the spin accumulation effect. 



H. EIGENFUNCTIONS AND EIGENVALUES OF THE ECCENTRIC QR WITH RASHBA SO 
INTERACTION 

We start by considering a 2D electron gas in the xy plane confined to an annular region delim- 
ited by two circles of radii a and b whose centers O and O' are separated by a distance d < b — a, 
that defines the eccentricity of the QR. A magnetic flux $ threads the structure (see FigQ])- 

The single particle Hamiltonian describing an electron of effective mass m* subjected to the 
Rashba spin orbit (RSO) coupling reads 

H = & + v + ><^->< (1) 

where a is the strength of the Rashba spin orbit (RSO) coupling and the Pauli matrices a are 
defined as standard. Let the polar coordinates be (p, 9) and (p', 9') with respect to O and O' 
respectively. The confining potential defining the eccentric QR is 



. for a < p < p ext = \b 2 + d 2 -2bd cos 9\/l — (v sin l 
V(p,9) = { H H V V Kb (2) 

oo otherwise . 

The vector potential which is introduced in the Hamiltonian via the substitution, p = Kk = 
—iKV — | A, is written in the axial gauge as A = (§/2np)9 . Using , b p = cos9a x + sm9a y and 
b Q = — sin 9a x + cos 9a y we can rewrite the Hamiltonian as 



H -. « r 



2m* 



-d p (pd p ) - ^{idg + v) 
p p 2 



+iacrgd p - -b p (id e + vj , (3) 



P 

where v = ^- is the magnetic flux in units of the flux quantum $ = hc/e. 
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FIG. 1: Schematic of the annular cavity with RSO threaded by a finite flux <£, where a and b are the internal 
and external radii respectively, and d the eccentricity. p(p'),6(6 ) are the polar coordinates referred to the 
origin 0(0 ). 

In what follows it will be useful to define the dimensionless coordinates £ = p/b, £' = p'/i>, 
the aspect ratio X — b/a and 

ma. 



h 2 E 

Due to the RSO, the bulk spectrum has two branches 17 

e= (fo/fc) 2 ± /5 (6fc) 



/; 2 



(4) 



(5) 



and for a given value of e there are two non-trivial solutions for the momentum k that we denote 
k + and k~ respectively. 

For the concentric ring (d = 0), the total angular momentum J z = l z + \ha z is a constant of 
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motion and j = m + \ is a good quantum number even for finite RSO coupling. In Ref. 
complete set of eigenfunctions (eigenspinors labeled by the quantum number j = m + 1/2) and 
eigenvalues have been found for the concentric QR with RSO. 

For a finite d > 0, the angular momentum is no longer conserved. In this case, we expand the 
solution of the Helmholtz equation associated to the hamiltonian Eq.© in a basis of eigenfunctions 
of the total angular momentum, expressed in polar coordinates (p, 9) referred to the origin O (see 
Fig.©, 



( 



\ 



y c 



3 i(m+l)0 I 



(6) 



<Pm-v+l (p) , J 



where ip m _ v (p) and <£ m _ u+ i (p) are linear combinations of Bessels functions of first and second 
kind J m - u (k + ~ p), J m -v+i{k + ~ p) and Y m ^ v (k + ~ p), Y m ^ v+1 (k + '~ p) that satisfy the boundary 
condition at the inner circle ip(p — a, 9) — 0; V0, i.e: 



^ m _ v (p) = A m (F' m _ u J m _ u {k + p) + F m ^Y m ^{k + p) + J m _ v {k-p)) (7) 
+ B m (G' m _ u J m ^{k + p) + G m ^Y m _ u {k + p) + Y m _ u {k-p)) 

VL-V+I (p) = An [F'm-v Jm-v+1 (k + p) + F m _ l/ Y m _ l/+l (k + p) - J m - U+1 (k~ p)) 

+ B m (G' m _ u J m ~ u+ i(k + p) + G m - u Y m ^ u+ i(k + p) — Y m ^ u+ i(k p)) . 

with the coefficients F m _ u , F' m _ v , G m _ u and G' m _ u given by : 

F m - U = D [Jm-v+i(k a)J m _ u ( y k + a) — J m -v(k a)J m — v +i(k + a)) 

F' m _ v = D (j m _ v+1 (k~a)Y m _ v (k + a) - J m _ v (k~a)Y m _ v+1 {k + a)) 

G m -u = D (Y m - u +i(k a)J m ^ u (k + a) — Y m _ v {k a)Y m — u +i(k + a)) 

G' m _ v = D (Y m - v+ i(k + a)Y m ^ u (k~ a) — Y m ^ v (k + a)Y m - u+ i(k~ a)) , 

D = (J m - u {k + a)Y m _ u+1 (k + a) - J m ^ u+1 (k + a)Y m _ u (k + a)) 1 . (8) 

The boundary condition ip(p ex t,9) = , V0 leads to a very complicated set of equation for 
determining the eigenvalues. One way to tackle this problem is to start by expanding the function 
ip(p, 9) into a Fourier series in 9'; 



^(p,9) 



jti— — oo Z— im=— oo ^mn^ 




where the matrices Q^ mn and Q}^ n are defined as : 

QL = - r Vl- V (p) e im0 e- ine 'd9' (9) 



The fully calculations of the above integrals are rather cumbersome since p and 6> are functions 
of p', 9'. To evaluate Eqs.© at the external boundary p = p ext , i.e p' = b, we employ addition 
theorems of Bessel functions 28 and solve the integrals analytically by means of a perturbative 
expansion to order 0(2) in the parameter v 5 < 1 (5 = d/b). After a lengthly derivation, we arrive 
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to the following set of equations : 



EE4" = 



rt 



fit 



(10) 



ft 



III 



where it is understood that — oo < n, m < oo being 



Vn 



mn 



mn 



(k + , J) + Fj 



Z mn {k + , y) + Z mn (k , </)) 



+ 5 m (G m _ u Z mn {k + , J) + G m ^ v Z mn {k + ,Y) + Z mn {k , F)) 
Q l mn = A m (F^W mn (k + ,J)+F m ^W mn (k + 1 Y)~W mn (k-,J)) 

+ fi m (G^W mn (fc + ,J) + W' m _ ! ,Z mn (A ; + ,r)-^ m „(A;-,r)) , (11) 



Z mn (k*,C) = a 1 C n - v (k%)J n _ m (k*d) - a 2 C n - v+1 (k*b)J n _ m+1 (k*d) 

+ a 2 C n - v - 1 (k*b)J n - m - 1 (k*d) + a 3 C n - u+2 (k*b) J n ^ m+2 (k*d) 
+ a A C n „ v _ 2 {k*b)J n ^ m _ 2 (k*d)) 

W mn (k*,C) = aiC n ^ u (k*b)J n ^ m ^i(k*d) - a 2 C n ^ u+1 (k*b)J n _ m (k*d) 

+ a 2 C n - u -i(k*b)J n - m - 2 (k*d) + a 3 C n - u+2 (k*b)J n „ m+1 (k*d) 

+ a A C n ^ 2 {k*b)J n ^{k*d)) . (12) 



In the last set of equation we denote k* = k + , k and the coefficients eij, z = 1, 2, 3, 4 are ai = 



2 - ^f-, a 2 = i/<J, a 3 = ^ (| - 1) and a 4 = (| + l). C = J, F is a Bessel function of first 



or second kind whenever it corresponds. 

For a given set of parameters /3, A and 5, we solve Eqs. (fT0l) - (fl~2~)) for different values of the 
flux v. The number of equations (defined by the maximum values of m, n) is increased until 
convergence in the solution is reached. The (dimensionless) energies e(/c + , k~) are found by the 
bisection method with a precision ~ 1CT 10 . We denote by the associated eigenspinor built with 
the expansion coefficients (An? B m ) that solve Eq.dTQ]) for a given e^. 

In order to fix numerical estimates for the parameters we consider characteristic values ex- 
tracted from experiments. Rings with external radius b ~ 400 — 500nm and an aspect ratio 
A ~ 2 have been recently employed as devices 7 . Typical values for the Fermi wavelength are 
X F ~ 40 — 50nm that give k F ~ O.lnm -1 . For an effective mass m* ~ 0.042m e , a Rashba cou- 



with the functions 



Z mn (k*, C) and W mn (k*, C) defined as: 
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pling constant a = 8meV ■ nm and b ~ AOOnm we obtain /3 ~ 4. These parameters characterize 
the sample S studied in the present work. 

In the case of concentric QRs (d = 0), the total angular momentum J is a constant of motion but 
the SO interaction breaks the degeneracy between states differing in one unit of j. The degeneracy 
between states with opposite values of j is removed by the finite magnetic flux v, being the charge 
persistent currents the signature of this broken time reversal symmetry 29 . The energy spectrum 
exhibits a structure of level crossings as a function of v, that are a signature of the Aharonov- 
Casher phase generated by the RSO coupling 16 . 

For d > the total angular momentum is no longer conserved due to the lost of rotational 
symmetry induced by the finite eccentricity. As a consequence, the quantum spectrum exhibits 
avoided crossings that are the manifestation of the dynamical tunneling between formerly degen- 
erated states 25 . 

We display in Figf2]the energy spectrum as a function of the flux v for d = 0.05 (black filled 
circles). This value of d is chosen to emphasize that a rather small eccentricity (~ 20nm for the 
sample S) is enough to generate qualitative changes in the quantum spectrum and its associated 
eigenfunctions, as we will show below. For comparison we show the spectrum for d = (turquoise 
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empty circles), already studied in Ref. 

The left panel of Figf2]shows the lowest eigenvalues of the QR S. Only the first transverse mode 
is active in this region and the spectrum shares some characteristics with the one of a ID ring (due 
to the symmetry respect to v = 0.5 the spectrum is shown for < v < 0.5). However, for 
d = 0.05 the lower levels are almost flat as a function of the flux, with associated eigenfunctions 
that are highly localized in one side of the QR. They are formed mainly by contributions from 
low angular momentum eigenstates of the concentric ring, being these states highly affected by 
the finite eccentricity. In Figj3]a) we plot a contour plot of the probability amplitude of the ground 
state l^ej 2 at v = 0.2, which it is mostly concentrated around 9 = ir. Notice that the finite 
eccentricity produces a similar effect than an in plane electric field, although it does not depend 
on the charge of the particle 26 . 

The other type of levels which lay higher in the spectrum have a finite slope det / dv (current 
carrying states) and resemble the ones of the unperturbed ring with higher values of the total angu- 
lar momentum j. However, due to the finite eccentricity, energy splittings and avoided crossings 
appear at v > 0. A blow up image of one of these avoided crossings is displayed in Figf2]in the 
box label a). Regarding the associated eigenfunctions, they are mainly generated by combination 
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FIG. 2: Color on-line. Dimensionless energies e as a function of the magnetic flux v for the annular cavity 
S with RSO coupling defined in the text. The spectrum for the concentric ring d = is plotted by turquoise 
open circles and for the eccentric QR, d = 0.05, by black dotted lines. Left panel: Lowest eigenvalues. 
Central panel: Intermediate energy eigenvalues displaying avoided crossings for d = 0.05. The regions 
inside the red boxes a, b, c are enlarged in the right panels to show details of the avoided crossings. 

of whispering gallery modes of the concentric ring 25 and therefore are not much affected by the fi- 
nite eccentricity. In Figj3]b) we display for v = 0.2 a contour plot of the probability amplitude for 
one of these states (e = 49.11) which is quite homogeneously distributed in the angular direction. 

When higher transverse channels are active, the spectrum for the eccentric QR is more involved. 
It displays additional avoided crossings generated by the mixing of levels that for d — were 
degenerated belonging to different transverse channels. Two detailed images of these avoided 
crossings are displayed in the right panel of Figj2]in boxes b and c. 

Again, localized (fiat levels) and current carrying states coexist in this region of the spectrum. 
As an example, in Figj3]c) we plot for v = 0.2 a contour plot of the probability amplitude for one 
of the angular localized states with e = 164.97. In displays maxima and minima along the radial 
direction which are the fingerprint of basis states belonging to the second transverse channel of 
the concentric QR. 

Taking into account that a given eigenspinor ip e . is an expansion in a basis of functions of 
total angular momentum, its probability amplitude has an angular profile which originates in the 

interference terms oc £ mm , (v>L-v ( pL>- v + ( pL-v+i ( Pin>-v+i) cos ( m ~ m ') d ( see E( l©- As we 
have shown in FigJH the angular profile depends on the specific eigenstate. This effect will also 
influence the spin accumulation describe in the next section. 
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FIG. 3: Contour plots of the probability amplitudes in the sample S for v = 0.2 and an eccentricity 
d = 0.05. a) ground state with e = 26.96, b) eigenstate with e = 49.11 and c) eigenstate with e = 164.97. 
The values of e are given at v = 0. 

III. SPIN ACCUMULATION EFFECT 

As a result of the RSO interaction, the spin projection is not a good quantum number. Given an 
eigenspinor tp €i whose general form is given by the Eqj6]the mean value of the z-projection of the 
spin density is proportional to 

M 

<r,(Z,8) = ^Ifz^ = Yl (l^-^l 2 - \V l m-u + l\ 2 ) + 

m=—M 

M M 

+ 2 E H (vl-^L,-„-vLv+ivL>-„ +1 )cos((m-m')e) , (13) 

m=—M m'=m+l 

where we have used the fact that the radial part of the up and down components of an eigenspinor 
are real functions. Thus the products (p]^_ v (p^, _ u are even under the interchange m m'. The 
sum is for — M < m < M, where M is given by the cut off in the expansion in the basis of total 
angular momentum eigenspinors. 

In Refllq it has been shown that in a multichannel concentric QR with RSO interaction, a 
spin accumulation effect (SAE) develops for finite values of the magnetic flux v. The SAE is the 
tendency of the total spin density for N particles, S z = J2n a z, to be different from zero, positive 
on one border of the sample and negative on the other one^. 

In the concentric geometry and for v = 0, states with opposite value of j have opposite values 
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FIG. 4: Left panelxontour plot of the spin density a z (£, 9) for the three eigenspinors displayed in FigfJ] 
The anisotropic pattern in the angular direction is clearly observed for a) and c), where for panel b) the 
angular profile is almost uniform. The right panels show a z vs £ for two specific angular directions selected 
in order to stress the effect. In all the cases is v = 0.2. 

of a z 7^ and then for even number of particles iV = 2p, is T, z = 0. A finite v ^ breaks the 
spatial symmetry between single particle states with opposite value of j. This can be understood 
in terms of the effective orbital index of the Bessel functions m e // = m — u, which turns to be 
m — v and — m — 1 — v for ji = m + 1/2 and = —ji, respectively. Thus the modulus of 
m e ff decreases for ji and increases for j 2 . As a consequence a z is "pushed" toward the internal 
(external) boundary of the sample for the state label by j\ (j'2). This is the main ingredient that 
contributes to generate the SAE in the concentric QRs ! -,/.e a total spin polarization Y> z that has 
opposite signs on each border of the sample S. However neither a z nor Y> z depend on the angle 9. 

In the eccentric QR, the SAE develops an anisotropic angular profile due to the interference 
(last) term in Eq{T3l The particular profile depends on the specific eigenspinor, exhibiting different 
characteristics for localized or current carrying states. As an example in FigjH we plot the contours 
plots of (T z (£, 9) calculated for the three eigenspinors analyzed in SecHO The anisotropic angular 
pattern is clearly observed for the two localized states a) and c), where it is almost uniform for the 
current carrying state b). 
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As the particle number iV is increased to relevant experimental values, transverse channels are 
activated and the description of the effect becomes more involved. Nevertheless, mainly due to the 
existence of localized states along all energy ranges, the anisotropy in the SAE is a generic feature 
of the eccentric QR pierced by a magnetic flux. 

As an illustration, in the left panels of FigJ5] we show the contour plots of the total spin density 
H z (£, 9) = °~z for iV = 100 particles and for two different values of the magnetic flux v = 
0.2,0.45. 

The sensitivity of the SAE to the value of the magnetic flux v, is due to the qualitative changes 
that the eigenspinors experience at the avoided crossings in the flux-energy landscape. For this 
particular filling, the anisotropy in the SAE is more pronounced for v = 0.2 than for v = 0.45. 

Besides the accumulation effect, the total spin density per unit area M z , is different from zero 
and its value depends on 9. Let the differential area be dA = (b 2 — a 2 )d9 (A = 0.75-7T& 2 is the 
area of sample S). It is easy to verify that M z {9) = | J 5 9) £ <i£. Taking b = OAfxm and 
for N = 100 particles, M z ranges from M z (9 = 0) = -0.3/im" 2 to M z (9 = 100°) = 0.48/im" 2 
for v = 0.2. For v = 0.45, M z is generally smaller, but the differences are still important and 
maxima between M z (9 = 0) = 1.33 x 10"V™" 2 and M z (9 = 180°) = -0.046/xm -2 . Thus, 
M z experiences important changes and even reverses it sign on distances of the order 1/im. This 
is scketched in the right panels of Figf51 where we show E 2 vs £ for the two angular directions 9 
especially selected to stress the anistropy. 



IV. DISCUSSION AND SUMMARY 



We have shown that a finite magnetic flux in an eccentric multichannel QR with SO interac- 
tion induces a spin accumulation effect (SAE) that exhibits an anisotropic profile in the angular 
direction. The effect, which is reminiscent to the drift of the spin polarization induced by the ap- 
plication of an in plane electric field, remains appreciable for filling numbers (N ~ 100) which 
are experimentally realized in actual semiconductor QRs. 

The spatial structure of the SAE which can be ~ 1/im in the sample S (see Fig©, is not far 
from the sensitivity of detection methods based on scanning Kerr microscopy, which has been 
recently employed to image spin polarization in semiconductor channels!-^. 

In addition, the total spin density per unit area M z , reaches at some angular directions, values 
~ 0.5/im~ 2 . These values are comparable to the spin densities measured so far in the Spin Hall 
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FIG. 5: Left panels: Contour plot of the total spin density 6) for N = 100 particles. Right panels: 
S z vs £ for the two angular directions 9 selected to stress the anistropy. Lower (upper) panels correspond 
to v = 0.2 (y = 0.45). 

Effect in semiconductor channels 18 . 

In the presence of an external magnetic field perturbation, the spin magnetization should be 
proportional to M z . With the help of new experimental techniques based on microresonators— or 
scanning magnetic force microscopes that can reach resolutions of ~ 10rm?3, it should be feasible 
to sense differences in the values of M z bewteen sample regions separated lfim. 

We believe that our results could help to manipulate and functionalize the spin polarization in 
QRs with SO interaction, without using applied electric fields or external currents. 
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